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Abstract 

Gessel's walks are the planar walks that move within the positive quadrant by 
unit steps in any of the following directions: West, North-East, East and South- West. 
In this paper, we find an explicit expression for the trivariate generating function 
counting the Gessel's walks with k ^ steps, which start at (0, 0) and end at a given 
point (i,j) G 
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1 Introduction 

The enumeration of lattice walks is a classical problem in combinatorics and this article is 
about the special case of Gessel's walks. These are the planar walks that move within the 
positive quadrant by unit steps in any of the following directions: West, North-East, 
East and South- West. To be more precise, let the walk be of length k and end at 
£ if i, j > then the next step can either be at (i — 1, j), (i + 1, j + 1), (i + 1, j) 
or (i — — 1); if i > and j = it can be at (i — 1, 0), (i + 1, 1) or (i + 1, 0); if i = and 
j ^ it can be at (1, j + 1) or (1, j). This is illustrated on Figure [T] below. 
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Figure 1: Steps in Gessel's walks 
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For e 7L\ and k ^ 0, let 

q(i,j,k) = #{ Gessel's walks of length k starting at (0,0) and ending at (i, j)}. 

These walks are named after Ira Gessel, who conjectured in 2001 that for any k 0, 
g(0,0,2A;) = 16 fe [(5/6) fc (l/2) fc ]/[(2) fc (5/3) fc ], where (a) k = a{a + 1) • • • (a + k - 1). This 
was proven in 2009 by Kauers, Koutschan and Zeilberger [7]. 

Let Q(x,y,z) be the trivariate generating function counting Gessel's walks: 



Q{x,y,z) = ^2q(i,j,k)x'u-' : 



This series is entirely characterized by the generating functions Q(x, 0, z), Q(0,y,z) and 
Q(0,0, z) — which count the walks that end on the borders of the quadrant — via the 
following functional equation: 

xyz ( — I Vx + xy J Q(x,y,z) = zQ(x,0,z) + z(y + l)Q(0,y,z) - zQ(0,0, z) - xy. 

\x xy z J 

(1) 

This equation, a classical result [H [2], is valid a priori in the domain {(x,y,z) : |x| < 
l,jy| < l,|z| < 1/4}. 

In this article, with the help of complex analysis methods we obtain closed expressions 
for Q(x,0,z), Q(0,y,z) and Q(0, 0, z) from exploiting the functional equation ([T]). 

By using computer calculations, Bostan and Kauers [1] showed that Q(x, y, z) is 
algebraic and found the minimal polynomials of Q(x,0,z), Q(0,y,z) and Q(0, 0, z) — 
by "minimal polynomial" of an algebraic function F(u,z), we mean the unique monic 
polynomial with coefficients in C[u, z] dividing any polynomial which vanishes at F(u, z). 
Van Hoeij pQ then managed to compute Q(x, 0, z), Q(0, y, z) and Q(0, 0, z) as their roots. 

Furthermore, Gessel's model is one of the 2 8 models of walks in the quarter plane 
starting from the origin and allowing steps in the interior according to a given subset of 
{ — 1,0, l} 2 \ {(0,0)}. Bousquet-Melou and Mishna [2] provided a systematic analysis of 
the enumeration of these walks. After eliminating trivial models and those already solved 
[3] by reduction to walks in a half-plane (a simpler problem), 79 inherently different 
models remained. Following an idea of [4], they associated each of them with a group of 
birational transformations of C 2 — for details about this group see Subsection 15.11 of this 
article. The group is finite in 23 cases and infinite in the 56 others. Through functional 
equations analogous to ([T]), they found the generating functions for 22 of the models with 
an underlying finite group. The following property has been crucial to their analysis: in 
the (half-)orbit sum of the trivariate generating function, all terms, except for the one 
corresponding to the identity element of the group, have a positive part in x or in y equal 
to zero. This last property is not valid for Gessel's walks and consequently, Gessel's model 
is the only one with a finite group that is not solved in [2]. 

Our method of applying ([T]) is notably different and can be generalized to all 79 walks 
described above |10j . It heavily relies on the analytic approach developed by Fayolle, 
Iasnogorodski and Malyshev |4j and proceeds by reduction to boundary value problems 
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(BVPs) of Riemann-Carleman type. In the rest of the introduction we sketch this approach 
and explain how we adapt it to the enumeration of Gessel's walks. 

The authors of [3] find explicit expressions for the generating functions of the stationary 
distributions (^%j)ij^o of some ergodic random walks in the positive quadrant Z+, 
supposed to have four domains of spatial homogeneity: the interior : i,j > 0}, 

the real axis {{i,0) : i > 0}, the imaginary axis {(0,j) : j > 0} and the origin {(0,0)}. In 
the interior, the only possible non-zero jump probabilities correspond to the eight nearest 
neighbors. 

First, they reduce the problem to solving the following functional equation on {(x, y) : 
M < 1, [2/1 < 1}: 

K(x, y)U(x, y) = k(x, y)vr(x) + k(x, y)vr(y) + k (x, y)ir 0fi , (2) 

where the polynomials K(x, y), k(x, y), k(x, y) and ko(x, y) are known, while the functions 

n(x,y) = E» J >i 7r t,j a; '~V _1 J n ( x ) = J2i^iHo xl ~ l as wel1 as n(y) = E J >i 7r o,jy- J_1 are 
unknown but holomorphic in their unit disc; the constant 7ro 5 o is unknown as well. 

Second, they continue the functions tt(x) and 7?(y) meromorphically (with poles that 
can be identified) to the whole complex plane cut along some segments, see Chapter 3 of 

a- 

Third, they prove that both unknown meromorphic functions tt(x) and ir(y) satisfy 
boundary value conditions of Riemann-Carleman type, see (5.1.5) on page 95 of [4j. Using 
information on the poles of tt(x) and 7r(y), they reduce the problems to finding some new 
holomorphic functions that are solutions to BVPs of the same type, see pages 119-124, 
and particularly (5.4.10). If the index x of these BVPs is non-negative (which actually 
is the generic situation considered in [3]), their solutions are not unique but depend on 
X + 1 arbitrary constants (for the notion of the index see (5.2.7) on page 98 or (5.2.42) 
on page 108). Consequently, in Part 5.4, the authors of [3] reduce the BVPs above to 
finding holomorphic functions satisfying some new BVPs of Riemann-Carleman type with 
an index % = — 1. Finally, these last problems are uniquely solved by converting them into 
BVPs of Riemann-Hilbert type, see Theorem 5.2.8 on page 108 of |4j. The functions ir(x) 
and n(y) can be reconstructed from these solutions; the constant 7To,o is computed from 
the fact that Ejj^o^ii = 1- 

Compared to ([2]), our equation ([I]) seems somewhat more difficult to analyse, since it 
involves an additional parameter z. On the other hand, the unknowns zQ(x,0, z), z(y + 
l)Q(0,y, z) and zQ(0,0,z) of (TlJ have constant coefficients, unlike vr(x), tt(tj) and 7ro,o in 
(|2|). This fact implies (see Subsection [6|) that zQ(x,0,z) and z(y + 1)Q(0, y, z) can be 
continued to whole cut complex planes as holomorphic and not only meromorphic 
functions. It also entails that zQ(x, 0, z) and z{y + 1)Q(0, y, z) satisfy BVPs of Riemann- 
Carleman type with an index x = 0; whose solutions are unique, up to additive 
constants. Then, unlike |3J, we don't transform these problems anymore but we solve 
them directly. Their solutions uniquely determine the functions zQ(x,0,z) — zQ(0,0,z) 
and z{y + 1)Q(0, y, z) — zQ(0, 0, z). The quantity Q(0, 0, z) is then found easily, e.g. from 
(UJ) by making the substitution (x, y, z) = (0, —1, z) — which is such that the left-hand side 
of ([!]) vanishes — see also Remark [8j Finally, Q(x,y,z) is determined via ([T]). 
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2 Reduction to boundary value problems of Riemann- Carle- 
man type 

Assumption. In the sequel, we will suppose that z is fixed in ]0, l/4[. 

Before we state our main results, we must have a closer look at the kernel L(x, y, z) = 
xyz [l/x + 1/ (xy) + x + xy — 1/z] that appears in ([1]) and introduce some notations. The 
polynomial L(x, y, z) can be written as 

L{x, y, z) = a(y, z)x 2 + b(y, z)x + c(y, z) = a(x, z)y 2 + b(x, z)y + c(x, z), 

where a(y,z) = zy(y + 1), b(y,z) = —y, c(y,z) = z(y + 1) and a(x,z) = zx 2 , b(x,z) = 
zx 2 — x + z, c(x, z) = z. Define also 

d(y, z) = b(y, z) 2 - 4a(y, z)c(y, z), d(x, z) = b(x, z) 2 - 4a(x, z)c(x, z). 

We have L(x,y,z) = if and only if [b(y,z) + 2a(y,z)x] 2 = d(y,z) or equivalently 
[b(x, z) + 2a(x,z)y] 2 = d(x,z). In particular, the algebraic functions X(y,z) and Y(x,z) 
defined by L(X(y, z),y, z) = and L(x,Y(x, z), z) = respectively have two branches, 
namely 

X (y,z) = [-b{y,z) + d(y,z) l l 2 ]/[2a{y,z)\, X^z) = [-b(y, z) - d(y, z) 1 ' 2 ]/^, z)}, 
Y (x, z) = [-b(x, z) + d(x, z) 1/2 ]/[2a(x, z)], Y^x, z) = [-b(x, z) - d(x, z) l/2 ]/[2a(x, z)]. 

For any z G]0, l/4[, the polynomial d has one root equal to zero, say yi(z) = 0, as well 
as two real positive roots, that we denote by y2{z) = [1 — 8z 2 — (1 — 16z 2 ) 1//2 ]/[82: 2 ] and 
yz [z) = [1 - 8z 2 + (1 - 162f 2 ) x / 2 ]/[82f 2 ]. We have < y 2 (z) < 1 < y 3 (z). We also note 
Vi(z) = oo. The points yk(z), k £ {1, ... ,4} are the four branch points of the algebraic 
function X(y, z). 

Likewise, for all z E]0, l/4[, d has four real positive roots, that we denote by x\(z) = 
[l+2z-(l+4z) 1 /2]/[ 2z ], X2 ( z ) = [l-2z-(l-4z) 1 /2]/[ 2z ], X3 ( z ) = [l-2z+(l-4z) 1 /2]/[ 2z ] 

and 2:4(2) = [1 + 2z + (1 + 4z) 1 / 2 ]/[2z]. We have < Xi(z) < x 2 {z) < 1 < x 3 (z) < x A (z). 
The points Xk(z), k € {1,...,4} are the four branch points of the algebraic function 
Y(x,z). 

We now present some properties of the two branches of both X(y, z) and Y(x, z). 
Lemma 1. The following properties hold. 

(i) Xk(y,z), k G {0,1} are meromorphic on C \ ([yi(z), yi (z)] U [2/3(2), j/4 (z)]). On the 
latter domain, Xq has a simple zero at —1, no other zero and no pole; X\ has a 
simple pole at —1, no other pole and no zero. Moreover, both Xq and X\ become 
infinite at y\(z) = and zero at y^{z) = 00. 

(ii) For all y £ C, we have \Xo(y, z)\ < \Xi(y, z)\. 

(Hi) Yk(x,z), k G {0,1} are meromorphic on C \ ([xi(z),X2(z)] U [x 3 (z), 2:4(2)]). On the 
latter domain, Yq has a double zero at 00, no other zero and no pole; Y\ has a double 
pole at 0, no other pole and no zero. 
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(iv) For all x G C, we have \Yq(x, z)\ < \Y\(x, z)\. 

Neither X^{y^z), k G {0,1} is defined for y in the branch cuts [yi(z), 2/2(2)] and 
[2/3 (2) ,2/4 (2)]. However, the limits X^(y,z) defined by X£(y,z) = limX k (y,z) as y -> y 
from the upper side of the cuts and XZ(y, z) = limX^(y, z) as y — > y from the lower side 
of the cuts are well defined. These two quantities are each other's complex conjugate, since 
for y in the branch cuts, d(y, z) < 0. A similar remark holds for Yj~{x, z), k G {0, 1} and 
x in [xi(z),x 2 (z)\ or [x 3 (z), x 4 (z}]. Precisely, for y € [y 1 (z),y 2 (z)] and x G [xi(z), x 2 (z)], 
we have 



-6(y,z) =fi[-d(y,. 
2a(y,z) 



1 1/2 



y ± (x,z) 



-fe(x, z) =f i[ — d(x, 2) 
2a(x, z) 



1/2 



(3) 



Xf{y, z) = X([(y, z) and Y x (x, z) = Y^(x, z). Furthermore, the identities ([3]) are true for 
y G [2/3(2), 2/4(2)] and x G [23 (z), ^4(2)] respectively if one exchanges X^(y,z) in Xj^(y,z) 

and l"o ± ( x ' z ) m ^o F ( x ' z )- 

In order to state our main results, we also need the following lemma. 

Lemma 2. The curves X([yi(z), 2/2(2)], z) and Y([xi(z), x 2 (z)], z) satisfy the following 
properties. 

(i) They are symmetrical w.r.t. the real axis and not included in the unit disc. 

(ii) X ( [2/1 (z), 2/2(2)], z) contains 00 and K([xi(z), x%(z)], z) is closed. 

(Hi) They split the plane C in two connected components. We call <£X ( [y\ (z), 2/2(2)], z) 
and §fy([xi(z), ^2(2)], z) £/ie ones of the branchpoints x\(z) and yi(z) respectively. 
They are such that [xi(z), x 2 (z)] C &X([yi(z),y 2 (z)], z) C C \ [2:3(2), 2:4(2)] as well 
as [?/i (z), 2/2(2)] c3T(Mz),x 2 (z)],z) C C\ [2/3(2), 2/4(2)]. 



X([ yi (z),2/2(2)],2 



y([xi(z),x 2 (z)],z 



£1(2) x 2 (2) ^£3 (z)x 4 (z) 




Figure 2: The curves X ( [2/1 (z), 2/2(2)], z) and Y([xi(z), #2(2)], z) 

The proofs of Lemmas Q] and [2] are given in Part 5.3 of [I] for z = 1/4; they can be 
plainly extended up to z g]0, l/4[. 
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We shall soon see that zQ(x,0,z) and z(y + 1)Q(0, y, z) satisfy BVPs of Riemann- 
Carleman type. It will turn out that the underlying boundary conditions for zQ(x, 0, z) and 
z(y + l)Q(0,y,z) hold formally on the curves X([yi(z),y 2 (z)], z) and Y([x\(z), X2(z)], z) 
respectively. However, these curves are not included in the unit disc, see Lemma [2] above, 
therefore the functions zQ(x, 0, z) and z(y + 1)Q(0, y, z) are a priori not defined on them. 
For this reason we need, first of all, to continue the generating functions up to these curves. 
This is why we state the following theorem; its proof is postponed to Section [6J 

Theorem 3. The functions zQ(x,0,z) and z{y + 1)Q(0, y, z) can be continued as 
holomorphic functions from their unit disc up to C \ [2:3(2;), x^(z)] and C \ [2/3(2), y±{z)\ 
respectively. Furthermore, for any y G C \ ([7/1(2;), 7/2(2)] U [2/3(2:), 7/4 (z)]), 

zQ(X (y, z), 0, z) + z(y + 1)Q(0, y, z) - zQ(0, 0, z) - X (y, z)y = (4) 

and for any x G C \ ([xi(z) 1 X2(z)} U [#3(2), X4 (z)]), 

zQ{x, 0, z) + z(Y (x, z) + 1)Q(0, Y (x, z), z) - zQ{0, 0, z) - xY (x, z) = 0. (5) 

It is immediate from Theorem [3] that for any z g]0, l/4[, the function Q(0, y, z) can be 
continued as a holomorphic function from the unit disc up to C\ [2/3(2), 2/4(2)] as well: the 
point y = — 1 cannot be a pole of Q(0, y, z) since the series X^>o fe^o 9(0> i> k) zk converges. 

Now, we derive the above-mentioned boundary conditions satisfied by the functions 
zQ(x, 0, 2) and z(y + 1)Q(0, y, z) on the curves X( [2/1(2), y 2 (2)], z) &ndY([xi(z),x 2 (z)],z) 
respectively. 

Lemma 4. zQ(x,0,z) and z(y + l)Q(0,y, z) belong to the class of the functions holomor- 
phic in ^X([yi(z),y 2 (z)], z) and &Y( [xi(z), ^2(2)] , 2) respectively, continuous up to the 
boundary of the latter sets and satisfying the boundary conditions 

z[Q(t,0,z)-Q(t,0,z)] = tY (t,z)-tY (t,z), Vt G X{[ Vl {z), y 2 (z)}, z), 

z[(t + l)Q(0,t, z) - (t + l)Q(0,t, z)} = X (t,z)t-X (t,z)t, Vt G Y([x 1 (z),x 2 (z)],z). 

(6) 

Proof. Due to Lemma [2] and Theorem [31 it remains to prove Equation ([6]) above. Let 
y G [yi(z) , y 2 (z)] and let y + and y~ be close to y, such that y + is in the upper half-plane 
and y~ in the lower half-plane. Then we have ([!]) for both y + and y~. If now y + — > y and 
y"~ ->■ y, then we obtain X (y + , 2;) -> 2) and Xo(y _ , 2) ->■ X^(y, z) = Xf(y, z). So 

we have proved that for any y G [1/1(2;), 1/2(21)], 

2Q(X+(y, 2), 0, 2) + z(y + 1)Q(0, y, 2) - zQ(0, 0, 2) - X+(y, z)y = 0, (7) 
zQ(X+(y, x),0, z) + z(y + 1)Q(0, y, 2) - zQ(0, 0, 2) - X+(t/, z)y = 0. (8) 

Subtracting (jSJ) from ([7]) gives that for 7/ € [2/1(2), 7/2 (z)], 

z[Q(X+(y, z),0, z) - Q{X+(y, z),0, z)] = X+(y, z)y - X+(y, z)y. 

Then, using the fact that for k G {0, 1}, y G [2/1(2), 7/2(2)] and 2 g]0, l/4[, Y (X^(y, z), z) = 
y — which can be proved by elementary considerations starting from Lemma [l] — we get the 
first part of ©. Likewise, we could prove the second part of ©. □ 
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3 Results 



Problems as in Lemma H] are usually called BVPs of Riemann-Carleman type, see e.g. 
Part 5.2.5 of A standard way to solve them consists in converting them into BVPs 
of Riemann-Hilbert type (i.e. with boundary conditions on segments) by using conformal 
gluing functions (CGFs), as in Equation (17.4') on page 130 of [5j. 

Definition 5. Let ?cCU {oo} be an open and simply connected set, symmetrical w.r.t. 
the real axis and different from 0, C and C U {oo}. A function w is said to be a CGF for 
set if: 

(i) w is meromorphic in '; 

(ii) w establishes a conformal mapping of& onto the complex plane cut along some arc; 
(Hi) for all t in the boundary of&, w(t) = w(t). 

Let w{t,z) andw(M) be CGFs for @X([yi(z),y 2 (z)],z) and &Y([xi(z), x 2 (z)}, z). The 
existence (but no explicit expression) of w and w follows from general results on conformal 
gluing, see e.g. Part 17.5 in [SJ. 

Transforming the BVPs of Riemann-Carleman type into BVPs of Riemann-Hilbert 
type thanks to w and w, solving the latter and working out the solutions, we will prove 
the following. 

Theorem 6. For z g]0, l/4[ and x € C \ [xs(z), Xi(z)], 
z[Q(x,Q,z) -Q(0,0,*)] = 



t , 1 
xY (x,z) -\ — 

7T 



»W t[-d(t,z)] 

xi(z) 



1/2 



d t w(t, z) 



d t w(t,z) 



2a(t,z) 

where w is a CGF for the set @X([yi(z),y2(z)],z). 

For z £]0, l/4[ and y £ C \ \y 3 (z),Vi(z)], 
z[(y + l)Q(0,y,z)-Q(0,0,z)] = 

~ -, 1 /O 

d t w(t,z) 



w(t, z) — w(x, z) w(t, z) — w(0, z) 



dt. 



1 rwW t[-d(t,z)] 1/2 

X (y,z)y+- 

l yi (z) 2a{t,z) 



IT 



d t w(t,z) 



w(t, z) — w(y, z) w(t, z) — w(0, z) 

where w is a CGF for the set &Y([xi(z),X2(z)], z). 
For z G]0,l/4[, 



dt, 



Q(Q,0,z) 



1 /■»(*) t[-d{t,z)] 
yi (z) 2a{t,z) 



7T 



1/2 



d t w(t,z) 



d t w(t,z) 



w(t, z) — w(— 1, z) w(t, z) — w(0, z) 



dt, 



where w is a CGF for the set %?Y([xi(z),X2(z)], z). 

The function Q(x, y, z) has then the explicit expression obtained by using the ones of 
Q(x,0,z), Q(0,y,z) and Q(0, 0, z) in |7]). 
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All functions in the integrands above are explicit, except for the CGFs w and w. In [3], 
suitable CGFs are computed implicitly by means of the inverse of some known function, 
see Equations (|25p and (|27[) in Section for the details. Starting from this representation, 
we are able to make explicit these functions for the case of Gessel's walks. In order to 
state the result, we need to define 



G 2 (z) = (4/27)(l + 224z 2 + 256z 4 ), 



(9) 



G 3 (z) = (8/729)(l + 16z 2 )(l-24z + 16z 2 )(l + 24z + 16z 2 ). 

We define also K(z) as the unique real positive solution to 

K 4 - G 2 (z)K 2 /2 - G 3 (z)K - G 2 (z) 2 /A8 = 0. (10) 

With the notations r k (z) = [G 2 {z)-exp(2knr/3){G 2 (z) 3 -27G 3 {z) 2 ) 1 ^}/3 for k € {0,1,2}, 
we have K(z) = [-ro(z) 1 / 2 + r x {zfl 2 + r 2 (z) 1/2 ]/2. We finally define 



1 - 24z + 16z 2 4(1 - 4z) 2 t 2 

r (t, Z) — 



3 z (t-x 2 (z))(t-l) 2 (t-x 3 (z))' 

l-24z + 16z 2 4(1 -4z) 2 t(t + l) 2 



(11) 



F(t,z) = + 



[(t-x 2 (z))(t-x 3 (z))} 2 



Theorem 7. A suitable CGF for the set @X([yi(z), y 2 (z)], z) is the unique function having 
a pole at x 2 (z) and solution to 

w 3 - w 2 [F(t,z) + 2K(z)] +w[2K(z)F(t,z) + K(z) 2 /3 + G 2 (z)/2] 

- [K(z) 2 F(t,z) + 19G 2 (z)K(z)/ 18 + G 3 (z) - 46K(z) 3 / 27} =0. 

(12) 

Likewise, a suitable CGF for the set&Y([xi(z),x 2 (z)],z) is the unique function having a 
pole at x 3 (z) and solution to the equation obtained from IW\) by replacing F by F , see 177]) . 

The function Q(x, y, z) being found in the domain {\x\ < 1, \y\ <l,z G]0, l/4[} thanks 
to Theorems [6] and [3 its coefficients Qij(z) = ^2 k>0 q{i,j, k)z k can be obtained most easily 
by the use of Cauchy formulas, for any real z G]0, l/4[. Since the radius of convergence 
of the complex series ^k^oli^ii^)^ ls n °t smaller than 1/4, the numbers q(i,j,k) of 
Gessel's walks can then be identified e.g. in terms of the limits of the successive derivatives 
as z — > 0+ (i.e. as z > goes to 0): 

Let us now outline two facts about the expressions of Q(x, 0, z), Q(0, y, z) and Q(0, 0, z) 
stated in Theorems [6] and [71 

Remark 8. The fact that ([TJ is valid at least on {|x| < 1, \y\ < 1, \z\ < 1/4} gives that for 
any triplet (x,y,z) G {\x\ < 1, \y\ < 1, \z\ < 1/4} such that L(x,y,z) = 0, the right-hand 
side of ([I]) equals zero, in such a way that 

z[Q(x,0,z) -Q{0,0,z)] + z[(y + l)Q(0,y, z) - Q(0,0, z)] + zQ{0, 0, z) - xy = 0. (13) 
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We deduce that 



zQ(0, 0, z) = -z [Q(x, 0, z) - Q(0, 0, z)] - z[{y + 1)Q(0, y, z) - Q(0, 0, 2)] + xy, (14) 



where the functions in square brackets in the right-hand side of (jl4|) are given thanks to 
Theorem [6l 

To get the expression of zQ(0, 0, z) given in Theorem El we have chosen to substitute 
(x, y, z) = (0, —1, z) in (|14p . which is suitable, since with Lemma[T]we have Xq(—1, z) = 0. 

Moreover, a consequence of Theorem [3] is that (113f) is valid not only on {L(x,y,z) = 
0} n {|x| < 1, \y\ < 1,2 e]0, l/4[} but in a much larger domain of the algebraic curve 
{L(x,y,z) = 0}. Namely, if (x,y,z) is such that z £]0, l/4[ and y = Yq(x,z) or x = 
Xo(y, z), then (|13|) is still valid. Substituting any triplet (x, y, z) lying in this domain into 
(Q3D yields zQ(0,0,z) as in {H). 

Remark 9. In Theorem z[Q(x,0,z) -Q(0,0,z)] and + 1)Q(0, y, -Q(0,0,z)] 
are written as the sums of two functions not holomorphic but singular near ^(z)] 
and [yi(z),y2(z)] respectively. The sums of these two singular functions are of course 
holomorphic near these segments, since the latter are included in the unit disc, according 
to Section [2) By an application of the residue theorem as in Section 4 of [8] , we could 
write both generating functions as functions manifestly holomorphic near these segments 
and having in fact their singularities near [2:3(2:), 2:4(2)] and [2/3(2), 2/4(2)] respectively. 

We conclude the discussion of Theorems El [6] and [7] with the following remark. 

Remark 10. With the analytical approach proposed in this article, it would be possible, 
without additional difficulty, to obtain explicitly the generating function of the number of 
walks of length k, starting at an arbitrary initial state (io,jo) and ending at Indeed, 
the only difference is that the product xy in (JT]) would then be replaced by x J ° +1 y JO+1 . 

The rest of the article is organized as follows. In Section HI we prove Theorem [H In 
Section \5\ we give the proof of Theorem [71 There, the general implicit representation of 
the CGFs inspired by [1] is given and developed for the case of Gessel's walks. The proof 
of Theorem [3] is postponed to the last Section [6l 

4 Proof of Theorem 6 

The proof is composed of three steps. 

Step 1. We solve the BVPs of Riemann-Carleman type with conditions (J6j) by transforming 
them into BVPs of Riemann-Hilbert type as in Parts 5.2.3-5.2.5 of [4] or in Part 17(.5) of 
[5]. The only notable difference from [4] is that the index of our problems is zero; this is 
why the solution zQ(x, 0, 2) is found in &X([yi(z),y2(z)],z) up to an additive function of 

2, as 




Vxe&X([ yi (z),y 2 (z)],z) 
(15) 
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where w is the CGF used for ^X([yi(z),y 2 (z)], z). Similarly, we could write an integral 
representation for z(y + 1)Q(0, y, z), up to some additive function of z. 

Step 2. We transform these representations into the integrals on real segments written in 
the statement of Theorem [6J This step is inspired by [8]. 
Let C(e, z) be any contour such that: 

(i) C(e, z) is connected and contains oo; 

(ii) C(e,z) C (&X([y 1 (z),y 2 (z)],z)UX([y 1 (z),y 2 (z)},z)) \ [x l (z),x 2 (z)}; 

(hi) lim e ^o C(e, z) = X([yt(z), y 2 (z)], z)l)S(z), with S(z) the segment [x 1 (z), X(y 2 (z), z)\ 
traversed from X(y 2 (z), z) to x\(z) along the lower edge of the slit and then back to 
X(y 2 (z),z) along the upper edge. 



X([ yi (z),y 2 (z)],z] 



xi(z) x 2 (z) ^xz{z) x 4 (z) 



C(e,z) 



Figure 3: The curve X([yi(z),y 2 (z)], z) and the new contour of integration C(e,z) 



Let £^C(e, z) be the connected component of C\C(e, z) which does not contain xs(z). Now 
we apply the residue theorem to the integrand of (I15p on the contour C(e, z). Thanks to 
Lemma Q] and Property (ii) above, t h > tlo(t,z) is holomorphic in 5^C(e, z). Likewise, 
by using Definition [5] and Property (ii) we reach the conclusion that dtw(t, z)/[w(t, z) — 
w(x, z)] is meromorphic on &C(e, z), with a single pole at t = x. For this reason, 

d t w(t,z) 



^- I tY (t,z) 



2iri 



C(e,z) 



w(t, z) — w(x, z) 



dt = xYq(x, z). 



(16) 



Then, letting e tend to 0, using Equations ()15p ~ (jl6p and Property (iii) of the contour, we 
derive that, up to an additive function of z, 



zQ(x, 0, z) = xY (x, z) 



1 

2ttl 



d tW (t, z) 

tY0{t, Z)— r r 

S(z) w{t,z) -w[x,z) 



dt. 



(17) 



Since for any x G &X([yi(z), y 2 (z)], z), the integrand of (fT7|) is, as a function of the variable 
t, holomorphic in ]x 2 (z), X(y 2 (z), z)[, we have 



I t y, (t ,.-) '] 

'S(z) W(t,z)-W{x, 



— dt= / ' [tY+(t,z)-tY-(t,z)] 

i z ) Jxi(z) 



d t w(t,z) 



w(t, z) — w(x, z) 



dt. 
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Using ©, we then immediately obtain the expression of z[Q(x, 0, z) — Q(0,0, z)] stated 
in Theorem [6] for x € &X([yi(z), 1/2(2)], z). Likewise, we could obtain the expression of 
z[(y + l)Q(0,y,z) - Q(0,0,z)} written in Theorem for y € WY([ Xl (z), x 2 (z)}, z). The 
formula for Q(0, 0,z) has already been proved in Remark [HI 

Step 3. In order to complete the proof of Theorem [6l we have to show that the integral 
representations of z[Q(x, 0, z) — Q(0, 0, z)\ and z[(y + 1)Q(0, y, z) — Q(0, 0, z)\ hold not only 
on 1fX([yi(z), 2/2(2)]) z) and ^Y([xi(z),x 2 (z)], z) but on the domains C\ [x^(z), x^(z)\ and 
C\ [2/3(2), 2/4(2)] respectively. 

It is clear that they can be continued up to C \ ([2:3(2;), 2:4(2)] U (w (w([xi(z),x 2 (z)], 



z)) \ [x 1 (z),x 2 (z)],z)) and C \ ([2/3(2), 2/4(2)] U (SH^Gs/i^), 2/2(2)], 2)) \ [2/1(2), 2/2(2)])) 



respectively. To conclude the proof of Theorem [6j it therefore suffices to show that 
w~ 1 (w([x 1 (z),x 2 (z)],z)) \ [ Xl (z),x 2 {z)} = and w~ l (w( [2/1 (2), 2/2(2)], 2)) \ [2/1(2), 2/2(2)] = 
0. This is the subject of Proposition [TTJ it is postponed to Section [5] because all necessary 



5 Study of the conformal gluing functions 

Notation. To be concise we drop, from now on, the dependence on z of all quantities. 

The main subject of Section [S] is to prove Theorem For this we shall define two 
functions, namely w in (j25j) and w in (|27p . which thanks to Part 5.5 of [3] are known to 
be suitable CGFs for the sets &X([yi,y 2 ]) and &Y([xi,x 2 ]) respectively. We shall then 
show that these functions satisfy the conclusions of Theorem [71 

5.1 Uniformization 

Let us begin Section [5] with studying a uniformization of the algebraic curve = {(x, y) £ 
(C U {oo}) 2 : L(x, y, z) = 0}, where L(x, y, z) = xyz [l/x + l/(xy) + x + xy — l/z] is the 
kernel appearing in ([T]). 

Proposition 11. For any z €]0, l/4[, J£ is a Riemann surface of genus one. 

Proof. We have shown in Section [2] that L(x, y,z) =0 if and only if [b(x)-\-2a(x)y] 2 = d(x). 
The Riemann surface of the square root of a polynomial which has four distinct roots of 
order one having genus one (see e.g. Part 4.9 in [BJ, particularly pages 162-163), the genus 
of J2? is also one. □ 

With Proposition [TTT it is immediate that Jzf is isomorphic to some torus. In other 
words, there exists a two-dimensional lattice O, such that Jz? is isomorphic to C/fi. Such 
a suitable lattice O, (in fact the only possible lattice, up to a homothetic transformation) 
is found explicitly in Parts 3.1 and 3.3 of [3J, namely uj\L + uj 2 Z, where 



facts about the functions w and w are proved there. 



□ 





(18) 
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We shall now give a uniformization of surface «£? , i.e. we shall make explicit x(ui) and 
y(uj), two functions elliptic w.r.t. lattice $7 and such that Jz? = {(x(uj),y(uj)),uj G C/Q}. 
By using the same arguments as in Part 3.3 of [1], we immediately see that we can take 



, n d'ixA 1 
x{u) =x 4 + —— — — — , y(uj) - 



p{uj)-d"{x i )/& yK ' 2a{x{u)) 



-b{x{u)) + d '^» 



2[p(u)-d"(x 4 )/Qf 



(19) 

where p is the Weierstrass elliptic function with periods cji,cj2- 

It is well-known (see e.g. (6.7.26) on page 159 in [9]) that p is characterized by its 
invariants g<i , g% through 

pV) 2 =4p(u;) 3 - 52 pH-<73- (20) 

Lemma 12. Invariants g2,93 of p are egua/ io: 

52 = (4/3) (l - 16z 2 + 16z 4 ), g 3 = -(8/27) (l - 8z 2 ) (l - 16z 2 - 8z 4 ) . (21) 

Proo/. We have that 4p(w) 3 - g 2 p{w) - g 3 = 4[p(uj)-p(uj 1 /2)}[p(uj)-p([uj 1 +uj2]/2)][p(uj)- 
p(w2/2)], see e.g. (6.7.16) on page 158 and (6.7.26) on page 159 in [9]. In particular, 
invariants 52,53 can be expressed in terms of the values of p at the half-periods. But it is 
clear by construction — and proved in Part 3.3 of p[| — that setting 

m = £M + fM t (22) 
6 t — Xi 

we have p(u>i/2) = /(X3), p{[oj\ + u>2]/2) = /(a^) and p(o>2/2) = f{x\). Lemma IT21 follows 
then from a direct calculation. □ 

For an upcoming use, we would like now to know the inverse images through the 
uniformization of the important cycles that are the branch cuts, that is to say x~ 1 ([xi, X2]), 
x _1 ([x3, X4]), y _1 ([yi, 2/2]) and y _1 ([y3, y^])- In this perspective, we introduce a new period, 
namely 

f Xl dx , , 

" 3 = Lri)F- (23) 

We will extensively use that LJ3 g]0, CJ2[ — this fact is proved in Lemma 3.3.3 on page 47 of 

m- 

Proposition 13. We have x~ 1 ([xi,X2]) = [0, uj\ [+U2/2 as well as x _1 ([x3, X4]) = [0, cji[, 
y~ l {[yi,y2\) = [0,wi[+[w 2 + uj 3 }/2 andy- l {\y^,y^) = [0, wi[+w 3 /2. 

Proposition 1131 follows from repeating the arguments given in Part 3.3 of [3]. It is 
illustrated in Figured! 

Let S(x,y) = 1/x + l/ (xy) + x + xy be the jump generating function of Gessel's walks. 
Consider the two birational transformations of (C U {oo}) 2 
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They satisfy ^ 2 = <J> 2 = id and Soty = So$ = S. Then, as in Part 2.4 of [4J, we 
define the group of the walk as group G generated by \& and <£. It is shown in [2] that for 
Gessel's walks, G is of order eight: in other words, inf{n > : ($ o ty) n = id} = 4. 

If (x,y) € (C U {oo}) 2 is such that L(x,y,z) = and if 8 is any element of G, then 
clearly L(9(x,y), z) = 0. This implies that group G can also be understood as a group of 
automorphisms of the algebraic curve It is then shown in (3.1.6) and (3.1.8) in Part 

3.1 of [4] that the automorphisms ^ and $ of become the automorphisms of C/fi 

ip(cu) = —uj, 4>{u) = —uj + u>3 (24) 

respectively. They are such that i/j 2 = (f) 2 = id, x o tp = x, y o ip = l/(x 2 y), x o <p = \/{xy) 
and y o <f> = y. 

A crucial fact is the following. 

Proposition 14. For all z s]0, l/4[, we have 0J3 = 3^2/4. 

Proof. Since the group generated by \& and $ is of order eight, the group generated by 
if} and cp is also of order eight, for any z e]0, l/4[, see Section 3 of [ID]. In other words, 
inf{n > : (0o^) ra = id} = 4. With (|24p . this immediately implies that 4co^ is some point 
of the lattice £1, contrary to W3, 2ojs and 3^3. But we already know that uj^ g]0, W2[, so 
that two possibilities remain: either W3 = W2/4 or W3 = 3^2/4. 

In addition, essentially because the covariance of Gessel's walks is positive, we can use 
the same arguments as in Section 4 of [8] and this way, we conclude that U3 has to be 
larger than W2/2, which finally yields Proposition 1141 □ 

5.2 Global properties of CGFs 

As said in Section EJ the existence of CGFs for the sets ^X([yi, 2/2]) and SfY"([xi, X2]) 
follows from general results on conformal gluing, see e.g. page 130 of Part 17.5 in [5]. 
Finding explicit expressions for CGFs is far more problematic: indeed, except for a few 
particular cases, like e.g. discs or ellipses, obtaining the expression of a CGF for a given 
set is, in general, quite a difficult task. 

But by using the same analysis as in Part 5.5 of [3], we obtain explicitly suitable CGFs 
for ( 3X{[yi,y2\) and @Y([xi, a^D- Before writing the expression of these CGFs, we recall 
that p, and thus also x with (I19p . take each value of CU{oo} twice on [0, W2[x [0, oj\/i[, but 
are one-to-one from the half-parallelogram (]0,W2/2[x [0, wi/&[)U[0, cji/2]U([0, wi/2] +1^2/2) 
onto C U {00} — indeed, see Corollary 3.10.8 in [6J and remember that p is even. In 
particular, on the latter domain, x admits an inverse function, that we denote by x^ 1 . 

Then, with Part 5.5.2.1 of [3], we state: 

w(t) = pi, 3 (z -1 (t) - [wi +W2]/2), (25) 

where is the Weierstrass elliptic function with periods 0)1,0)3, x^ 1 the inverse function 
of the first coordinate of the uniformization (|19f) and where wi, 6^2,^3 are defined in (|18l) 
and P}. With ([19]), we note that 

x- 1 (t) = p- 1 (f(t)), (26) 
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where / is defined in (|22|) . 

In Section 4 of [SJ, we have studied some properties of the function w defined in (|25p 
and we have shown that if u 3 > 002/2 (which is actually the case here, see Proposition I14p . 
then w is meromorphic on C \ [£3, £4] and has there a single pole, which is at x<i- 

Let us now notice that 

wit) = w(X (t)) (27) 

is a suitable CGF for the set &Y([xx,X2])- Indeed, as we verify below, the three items of 
Definition [5] are satisfied. 

Firstly, we immediately deduce from Lemma [TBI and from the inclusion &Y(\xi, X2]) C 



C\[ys, 2/4] stated in (iii) of Lemma[2]that w defined in (|27p is meromorphic on @Y([xi,X2\). 



Secondly, using that w is a CGF for @X([yi, yj[) as well as Property (i) of Lemma [T"5l 
we reach the conclusion that w establishes a conformal mapping of @Y([xx, X2]) onto the 
complex plane cut along some arc. 



Thirdly, once again with Property (i) of Lemma[[5j we get Xq(Yq(x)) = Xq(Yi(x)) = x 



for x € [xi,x 2 ]. As an immediate consequence, for t G ¥"([#1, x%\) we have w(t) = w(t). 

Lemma 15. The two following properties hold. 

(1) X :^Y([x 1 ,x 2 })\[y l ,y 2 }^^X([y ll y 2 })\[x 1 ,x 2 } and Y : &X([ yi ,y 2 ]) \ [x u x 2 ] -> 
< &Y([xx,x 2 \) \ [2/1,2/2] o,re conformal and inverse to one another. 

(ii) X (C) C C \ [x 3 , x 4 ] and Y (C) C C \ [y 3 ,y 4 ]. 

The proof of Lemma [T5l is done in Part 5.3 of [3] for z = 1/4; it can be generalized 
directly up to z e]0, l/4[. 



Lemma 16. The function w defined in (W)) is meromorphic on C \ [2/3,2/4] and has there 
a single pole, which is at x$. 

Proof. Xq being meromorphic on C \ ([2/1,2/2] U [2/3,2/4]) an d w on C \ [x 3 , x±], the function 
w defined in (j27| is a priori meromorphic on C\ ([2/1,2/2] U [2/3,2/4] UXq"^^,^])). But on 
the one hand, with (ii) of Lemma [TBI Xq ([a^,^]) = 0. And on the other hand, thanks 



to the gluing property of the CGF w, w satisfies w + (t) = w~(t) for t G [2/1,2/2], i.e. the 
limits of w(u) when u — > t £ [2/1,2/2] from the upper and lower sides of the cut are equal. 
w is thus also meromorphic in a neighborhood of [2/1,2/2], see e.g. Parts 5.2.3 and 5.2.4 of 
[3]. Finally, w is meromorphic on C \ [2/3,2/4]- 

Moreover, with ([271) and since w has on C \ [#3, X4] only one pole, which happens to be 
at x 2 , the only poles of w are at the points t where Xo(t) = x 2 . It is then easy to verify, 
by a direct calculation, that X3 is the only solution to the latter equation. □ 

We prove now the following proposition, which completes the proof of Theorem [6j 

Proposition 17. We have w~ 1 (w([x\, x 2 ]))\[xi, x 2 ] = and w~ l (w( [2/1, 2/2])) \ [2/1, 2/2] = 0- 

Proof. In order to prove the first identity, it is enough to show that for any fixed 
u € [sci,:c2]j the only solution in t of w(t) = w(u) is t = u. 
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If u £ [xi,X2], then w(u) 6 pi,3([— wi/2, wi/2]), see Proposition [131 and Equation ([25]) . 
Thus once again with (f25|) . the equation u)(i) = tu(u) can be interpreted as 

Pi, 3 (w) = Pi,s(T), (28) 

where w = x _1 (i) — [u>i + cl>2]/2 and T G [— wi/2, wi/2]. 

A priori, Equation (|28p admits the solutions u = ±T + fciCJi + A^a^, with fcj, k% € Z — 
see Corollary 3.10.8 in [6j, then remember that px s s is even and periodic w.r.t. uji,<jj$. But 
in our case, oj belongs to a restricted region, namely ] — 0J3, 0] x [— u\/{2i), oj%/(2i)]. Indeed, 
as already noted in this section, we have x _1 (CU{oo}) C [0, uj 2 /2] x [0, loi/i], so that thanks 
to Proposition [T31 we have x _1 (CU{oo}) - [uji + u 2 ]/2 C] -w 3 ,0] x [-uj 1 /(2i),oj 1 /(2i)]. 
In the latter restricted domain, the only solutions to (|28p are u = ±T. In particular, we 
get x _1 (i) - [oji + oj 2 }/2 = ±(x~ 1 (u) - [uji + u 2 ]/2), which yields x _1 (t) = (1 =F + 
u 2 ]/2±x- 1 (u). 

Then, taking the image of the previous equality trough p and using Equation (|26p . we 
obtain /(i) = p((l =F iM^i + ^2]/2 ± x~ 1 (n)). Since p is periodic w.r.t. coi,co 2 and even, 
we conclude that f(t) = f(u) and finally that t = u, since / is one-to-one, see ()22|) . 

The proof of the identity w (w([xi, x 2 ])) \ [x\,x 2 ] = is concluded. Similar reasoning 
and the use of Equation (|2T|) yield that ■w~ 1 (w([yi, y 2 ])) \ [2/1,2/2] = 0- D 

5.3 Proof of Theorem 7 

Our aim is to show that the function defined in (|25p is the only function having a pole at 
x 2 and solution to (|12p . 

Denote W4 = io 2 /4 and let pi 5 4 be the Weierstrass elliptic function with periods 011,014. 
We recall that p and px,3 are the Weierstrass elliptic functions with respective periods 
oji,u 2 and 011,013, where 013 = 3u> 2 /4 thanks to Proposition 1141 

To begin with, let us prove the following lemma. 

Lemma 18. Let p be the Weierstrass elliptic function with periods noted oj,oj and let n 
be a positive integer. Then the Weierstrass elliptic function with periods L0,Cb/n can be 
written in terms of p as follows: 

n-l 

+ [&( u + " P( fcc V n )] • ( 29 ) 

fc=i 

Proof. It is easy to verify that both the Weierstrass elliptic function having for periods 
£j,Lj/n and the function defined by (|29p satisfy the three properties hereafter: they are 
elliptic with periods u),u)/n; they have only one pole in the fundamental parallelogram 
oj[0, l[+(uj/n)[0, 1[, this pole is at and is of order two; they admit an expansion at to = 
equal to 1/lo 2 + 0(lu 2 ). Therefore, they must coincide, see e.g. Part 8.10 on pages 227-230 
in [9]. □ 

Now we notice that by applying the following addition formula (see (6.8.10) on page 
162 in [9j) 

p(u> +u) = -p(u) - p(u) + - 



p(u) - p(u) 



Vcj, oj, 



(30) 
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to the Weierstrass elliptic function p in ([2"9~]) and by then using the identity ([20]) . we can 
express the Weierstrass elliptic function with periods di,oi/n as a rational function of the 
Weierstrass elliptic function p with periods ui, 01. 

We shall apply this procedure in the proof of Lemmas [19] and [20] 

Lemma 19. We have 

(31) 

where p(ui2/2) = f(x\), p(w2/4) = (1 + 4z 2 )/3, p'(u;2/4) = —8z 2 and where p'(oi) can be 
expressed in terms of p(oi) and z with Equations [2U\) and [2l\) . 
Furthermore, 

P1,4(P _1 (/(*)) " [wi + W2]/2) = F(t), Vt € C, (32) 
where F is defined in [77]) and / in < fH] ). 

Proof. We have W4 = W2/4 by definition of 0J4. Then, with (|29p . we can write 

Pia(oi) = p(w) + p(a;+W2/2) + p(a;+a;2/4) + p(a; + 3u;2/4)-p(u;2/2)-p(w2/4)-p(3a;2/4). 

Using then addition formula ([30]) for p as well as the three equalities p(w2/4) = p(3w2/4), 
p'(o;2/4) = — p' (3w2/4) and p' (012/2) = — obtained from the facts that p(oi2/2 + 01) is 
even and p' (012/2 + 01) is odd, see (6.8.12) on page 162 in [S] — we get (|3ip . 
Using the formula below (see e.g. Exercise 8 on page 182 in [9 ) 

p(u 2 /A) = p(u 2 /2) + [(p(w2/2) - p(wi/2))(p(w2/2) - p([wi + wJ/2))] 1/2 (33) 

as well as p(wi/2) = f(xs), p([oi\ +W2]/2) = /(a^) and p(oi2/2) = f(x±), see the proof 
of LemmaCLl we immediately find p(w 2 /4) = (1 + 4z 2 )/3. With ([20]) and ([21]), we derive 
p'(cj2/4) 2 = 64z 4 . Since p is decreasing on ]0, W2/2[, see e.g. Part 6.11 on pages 166-167 
in [9], we have £>'(a>2/4) < and therefore ^(012/^) = — 8z 2 . 

Formula (]3T]) with the known values of p(u>2/2), p(w2/4), p'(w2/4) as well as with p'(ui) 
expressed in terms of p(ui) and z thanks to (]20p and (|2ip gives a representation of p\^(oi) 
as a rational function of p(oi). 

Evaluating this representation at 01 = p (f(t)) — [oi\ +W2]/2, once again using ([30]) for 
the function p together with (|20p and (I2ip for the derivatives as well as the explicit value of 
p([oi\ +W2]/2) given above, we get ([32]) . after a substantial but elementary calculation. □ 

Lemma 20. We have 

Pl,4(w) = -P1,3(W) + 7TT ^ ; 7 — 7^vi2 ~ 4 Pl,3(W3/3), Vw. (34) 

2[pi,3(w) - pi i3 (o;3/3)J 2 

Proof. Formulas ([29]) and ([30]) combined with the fact that 014 = 012/^ = 0)3/3, see 
Proposition HU easily lead to ([341 . □ 
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Equality ([20]) for pi )3 , written as 

p'1,3^) 2 = 4 Pi,3(^) 3 - 52,i,3pi,3(w) - 53,1,3, (35) 

allows us to express p[ 3 (w) 2 in terms of pi i3 (w) and invariants 52,1,3,53,1,3 associated with 
pi 3. The next lemma gives their expression in terms of z. 

Lemma 21. Invariants 52,1,3,53,1,3 0/ Pi,3 have the following explicit expressions: 

52,1,3 = 40p li3 (w 3 /3) 2 /3 - G 2 , 53,1,3 = -280p li3 (^ 3 /3) 3 /27 + 14p li3 (^ 3 /3)G 2 /9 + G 3 , 

where G 2 ,G 3 are defined in (OJ). 

Proof. The proof consists in expanding p\^{uj) at oj = in two different ways. 

Firstly, we use equality (f3T|) with p'(u) expressed in terms of p{oj) thanks to (f20|h with 
52,53 obtained in ([2T]) as well as with p(o>2/2), p(w2/4) and p'(a;2/4) found in Lemma [T9l 
Expanding this identity in a neighborhood of u = 0, we obtain: 

Pi, 4 M = 7J2 + [9G 2 /20]o; 2 - [27G 3 /28]u; 4 + 0(oj 6 ). (36) 

Secondly, we expand p 14 (a;) at w = using Equation (f34"|) with p' x 3 (w) and p' : 3 (w 3 /3) 
expressed thanks to f)35[) . After some calculation, we get: 

P M H = ^2 + [6pi, 3 (w 3 /3) 2 - 952,1,3/20] ^ 2 (37) 
+ [l0pi, 3 (^ 3 /3) 3 - 3pi, 3 (w 3 /3)52,i,3/2 - 2753,1,3/28] oj 4 + 0(lo 6 ). 

Lemma [2T1 follows then as we identify the expansions (|36p and (|3T|) . □ 

In the next lemma, we compute pi,3(ws/3). 

Lemma 22. W^e /iaue pi,3(w3/3) = K , where K is found as the only real positive solution 
to Equation U0\) . 

Proof. It is stated in Exercise 7 on page 182 in [9] that the quantity pi,3(w3/3) is the only 
real positive solution to K 4 — 52,1,3-K' 2 /2 — 5 3) i j3 -K" — g\ 1 3 /48 = 0. Replacing 5 2) i j3 and 
53,1,3 with their expression obtained in Lemma EH we conclude that pi j3 (o; 3 /3) is a root 
of' ' 

K 4 -G 2 K 2 /2-G 3 K -G 2 2 /48. (38) 

Let us now show that for any z G]0, l/4[, the polynomial (|38p has a unique real positive 
root. Since pi j3 (o; 3 /3) > — indeed, pi >3 is positive on [0,^3], see Part 6.11 on pages 166- 
167 in [9] — pi,3(w3/3) shall be characterized as the only real positive solution to (fTO)) , 

According to Lemma l23l it is now enoug h to verify that G 2 ^ and that G|-27G| > 0. 
The first fact is actually an immediate consequence of ([9]), while the second comes from 
the identity Gf, - 27G§ = (4 14 /3 6 )z 2 {z - l/4) 4 (z + 1/4) 4 , see also ©. □ 
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Lemma 23. For any real numbers G2,G^ such that G2 7^ and G 3 , — 27G| > 0, the 
polynomial l\38\) has a real negative root, a real positive root and two non-real complex 
conjugate roots. 

Proof. If some real numbers G2,G^ such that G2 7^ and G 3 , — 27G| 7^ are given, then 
there exists a lattice OfL + wZ and a Weierstrass elliptic function w.r.t. this lattice, say 
p, having for invariants G2,G%, see Corollary 6.5.8 on page 287 in [6]. By using once 
again Exercise 7 on page 182 in [9], we then come to the conclusion that p(w/3), p(tD/3), 
p{[uj — £D]/3), p([w + w]/3) are the four roots of the polynomial ([38]) . 

Now we prove that the inequality Gf, — 27G| > yields that one of p(u>/3), p(w/3) 
is negative while the other is positive, and that p([w — w]/3), p([w + &]/3) are complex 
conjugate to one another. 

If G\ - 27G§ > 0, then the periods qj,u) of p can be chosen such that ui > and 
uj/i > 0. Indeed, with pages 110-111 (particularly Theorem 3.6.12) of [6j, we conclude 
that 6b, Cb are either real and purely imaginary (if the polynomial 4x 3 — G2X — G3 has 
three real roots) or complex conjugate to one another (if 4x 3 — G2X — G3 has only one 
real root); in addition, it is well known that 4x 3 — G2X — G3 has three real roots if and 
only if Gl - T1G\ > and only one real root if and only if G\ — HG\ < 0. Moreover, 
on the parallelogram [0,o>] x [— Cb/(2i), tD/(2z)], p takes real values on the segments [0,cD], 
[— oj/2,u}/2], [0,u] ±w/2, [— cD/2, cD/2] +w/2, and [— a)/2,cD/2] + a) and on those segments 
only, see Part 3.16 on pages 109-115 in [6]. 



-cb/3 

-Lb/2 



, [^ + ^]/3. 

d>/3 


w/2 


1 — > 1 





Figure 4: The parallelogram [0, w] x [— u/(2i), to/(2i)] 

As a first consequence, p(d)/3) and p(tD/3) are real. Furthermore, on the segments 
[0,u] and [— cD/2,w/2], the values of p are positive and negative respectively, see Part 6.11 
on pages 166-167 in [9j. This is why p(w/3) > and p(w/3) < 0. 

As a second consequence, p([a) — cD]/3) and p([cD + £>]/3) are non-real. Since [u — w]/3 
and [u) + w]/3 are complex conjugate to one another, so are p([d) — cD]/3) and p([cu + cj]/3), 
see also Part 6.11 on pages 166-167 in [9]. Lemma [23] is proved. □ 

Remark 24. Lemma [23] is also valid for any real numbers G2,Gs such that G2 7^ and 
Gl - 27G| < 0. On the other hand, if Gf> - 27G^ = 0, then the polynomial ([38]) has only 
real roots. 
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Substituting in ([31)) equality (f35|) for p' 13 (u;) and p' x 3(^3/3), we express p\^{uj) in 
terms of pi,3(w), 52,1,3 , 53,1,3 and pi,3(w 3 /3). Applying then Lemmas I2T1 and [221 we get: 

Pi, 3 (w) 3 - Pi, 3 H 2 [pi,4(w) + 2#] + Pi, 3 H [2Jfpi, 4 (w) + ^ 2 /3 + G 2 /2] 

- [ET 2 pi,4(w) + 19G 2 K/18 + G 3 - 46is: 3 /27] = 0. 

In particular, evaluating this identity at oj = p~ 1 (f(t)) — [ui +U2}/2, using (p2|) and taking 
into account ([25)1. we obtain Equation ([12)1 for the CGF w; defined in ([25)h 

If i 7 , defined in (TlTI) . is infinite at some point to, then Equation (fT2l) yields [w(to) — 
JC] 2 = 0. Thus K is a double root of ([T2"D at to- I n addition, the double root K being 
non-zero by its definition via ([TDD and the product K 2 F(t ) + 19G 2 i^/18 + G 3 - A6K 3 /27 
of all roots of the polynomial (|12|) being infinite, the third root of (|l'2p must be infinite. 

Since F is infinite at x 2 , see ([TT)) . and since w has a pole at x 2 , see Subsection 15,21 
can thus be characterized as the unique solution to (|12p with a pole at x 2 — the two other 
solutions are finite at x 2 and equal to K. The part of Theorem [7] dealing with a CGF for 
the set &X([yi,y2]) is concluded. 

Now we prove the corresponding fact for &Y([xi, a? 2 ])- Since w is a solution to (|12h and 
since an easy calculation gives F(Xo(t)) = F(t), we conclude that w satisfies the equation 
obtained from (|12[) by replacing .F(i) by F(t). Furthermore, both F and w have a pole at 
X3: for i 7 , this is a consequence of (jlip and for w, this follows from Lemma[16l Using the 
same arguments as above for w, we then derive that w can be characterized as the only 
function having a pole at x 3 and solution to the equation obtained from (]12p by replacing 
F(t) by F(t). □ 

6 Holomorphic continuation of zQ(x, 0, z), z(y-\-l)Q(0, y, z) 

In this part, we shall prove Theorem[3l In other words, we shall show that zQ(x, 0, z) and 
z(y + 1)Q(Q, y, z) can be holomorphically continued from their unit disc up to C \ [x 3 , X4] 
and C \ [2/3 , 2/4] respectively. 

Proof of Theorem^ First of all, we lift the functions Q(x, 0, z) and Q(0,y,z) up to C/fi 
by setting q x (u)) = Q(x(u), 0, z) and ^(w) = Q(Q,y(u)),z). We recall that and y(u), 
the coordinates of the uniformization, are defined in (|19p . The functions and q y are 
a priori well defined on < 1}) and < 1}) respectively. Then, we use the 

following result, that we shall prove in a few lines. 

Theorem 25. q x [oS) and (y(u) + l)q y (uj), initially well defined on < 1}) and 

7/ 1 ({ 1 2/ 1 < 1}) respectively, can be holomorphically continued up to the whole parallelogram 
C/fi cut along [0,uj±[ and [0, uj\ [+w 3 /2 respectively. Moreover, these continuations satisfy 

q x (u) = <fo(V>(w)), Vw6C/f]\[0,u;i[, g v (w) = Vw € C/fi \ ([0, Wl [+w 3 /2), 

(39) 

and 

ZftcH + l)q y {uj) - zQ(0,0,z) - x{uj)y{uj) = 0, Vw G]3w 2 /8, w 2 [x [0, wi/«[. (40) 
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Finally, we set Q(x, 0, z) = q x (u>) if x{u) = x as well as Q(0, y, z) = q y (ui) if y(co) = y. 
Thanks to ([39]) and Proposition [131 these equalities define not ambiguously Q(x,0,z) and 
Q(0, y, z) on C\[xs, X4] and C\[?/3, 2/4] respectively, as holomorphic functions. Furthermore, 
the statements (j3J) and ([5]) are immediate consequences of PU|) . □ 



The proof of Theorem [3] is now reduced to that of Theorem [25] In order to carry out 
this proof, we need to find the location of the cycles a? -1 ({ |x| = 1}) and = 1}) on 

C/O. This is the subject of the following result, illustrated on Figure [5] below. 



X 4 



2/3 



2/1 



•7'2 



V2 



{\*\ < 1} 

{M<i} 



Figure 5: Location of the important cycles on the surface [0, w 2 [x [0, 

Proposition 26. We have = 1}) = ([0, Wi [+w 2 /4) U ([0, u>i [+3w 2 /4) as well as 

2T 1 ({M = 1}) = ([(W+u; 2 /8) U ([0,Wi[+5w2/8). 

Proof. The details are of course essentially the same for x and y, so that we are going to 
prove only the assertion dealing with x. 

First of all, we note that because of the equality x o ip = x, it is sufficient to prove that 
= 1}) PI ([0,u> 2 /2[x[0,u;i/i[) = [0, u;i [+w 2 /4 — the advantage of this is that p, and 
therefore also x, are one-to-one in the half-parallelogram [0, w 2 /2[x [0, qji/i[. 

The proof is then composed of three steps. 

Step 1. We prove that x(w 2 /4 + wi/2) = 1. For this, we recall that p(w 2 /4) = (l + 4z 2 )/3, 
p'(w 2 /4) = -8z 2 , p(wi/2) = /(x 3 ) with / defined in ([22]) and p'(wi/ 2 ) = 0, see LemmaCE] 
and its proof. With addition formula (j30p . we then immediately obtain the explicit value 
of p(w 2 /4 + uji/2). Finally, by using Equation (j!9p and after a simple calculation, we get 
x(w 2 /4 + wi/2) = 1. 

Step £ We show that = 1}) n ([0, w 2 /2[x [0, C [0, Wi[+w 2 /4. For this, let 

S [0,27r[ and consider the equation x{uj) = exp(i9). Thanks to ([T9]) and ([22]) . we obtain 
p{oj) = f(ex.p(i6)) and thus to = p _1 (/(exp(i#))). We can then use the explicit expression 
of the inverse function of p on [0, w 2 /2[x [0, uji/i[, see e.g. Part 6.12 on pages 167-172 in 
[9] , and we get 



w = x W + / 7173 1 ii/a = + o*/2 + - / — (41) 

Jf(l) W- - 92t - 93 V /A 2 JexoM) \d(x) \ L I Z 
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where d is defined in Section [2] and g 2 , 53 in Lemma [T2l Note that the second equality 
above comes from the first step and from exactly the same calculations as in Part 3.3 of 
[3]. Now we notice that d(x) = x A d(l/x). In particular, the change of variable x 1— > 1/x in 
the integral dx / [d(x)}^ yields dx/[d(x)\ V* = -J^diMi)] 1 / 2 . As 

a consequence, this integral belongs to iR. In conclusion, with (|41|) . we have shown that 
x-\{\x\ = 1}) n ([0,w 2 /2[x [0,wi/z[) C [0,wi[+W2/4. 

S'tep 5. We prove that the inclusion above has to be an equality. Indeed, if it was not 
the case, the curve = 1}) PI ([0, uj2/2[x [0, would be not closed, which is a 

manifest contradiction with the facts that {\x\ = 1} is closed and that x is meromorphic 
as well as one-to-one in the half-parallelogram [0, w 2 /2[x [0, uj\/i[. □ 

Proof of Theorem\25[ This proof is composed of two steps: at first, we shall define the 
continuations of q x and q y on the parallelogram C/fi cut along [0,wi[ and [0, wi[+a;3/2 
respectively; then, we shall verify that these continuations satisfy the conclusions of 
Theorem [ 



Step 1. We define the continuations of q x and q y . 

(i) We define q x (uS) on x _1 ({[x| < 1}) by Q(x(ui), 0, z) and q y (u) on < 1}) by 
Q(0,y(u), z). Note that as a consequence of Proposition \26\ we have < 
1}) = [w 2 /4,3w 2 /4] x [0,un/i[ and y-\{\y\ < 1}) = [5w 2 /8, 9w 2 /8] x [0,wi/*[. 

(ii) Motivated by ((TJ), on [3w 2 /4,a; 2 [x [0, Wi/i[c | < 1}), we set q x (u>) = -(y(u}) + 
l)q y (uj) + Q(0,0,z) +x(u)y(uj)/z and on }3ui 2 /8, 5w 2 /8] x [0,o;i/i[c a; -1 ({1^1 < 1}), 
we set (y(uj) + l)q y {uj) = —q x (u) + Q(0, 0, 2) + x(ui)y(uj)/z. 



(iii) On ]0,o; 2 /4] x [0,wi/z[, we define by q x ((f)(u))). Note that with Equation ((2 

we have </>(]0, w 2 /4] x [0,wi/z[) = [3w 2 /4, w 2 [x [0, wi/z[. On [w 2 /8, 3w 2 /8[x [0, 
we define q y (ui) by q y (ijj(ui)). By using ([2"3]) . we have ^([a; 2 /8, 3w 2 /8[x [0, 0Ji/t[) = 
]3w 2 /8,5w 2 /8] x [0,wi/t[. 

The functions g x and ^ are now well defined on the whole parallelogram C/Q cut along 
[0,o;i[ and [0, tJi[-kJ3/2 respectively. 

Step 2. We prove that the continuations of q x and q y defined in the first step satisfy the 
different assertions of Theorem [ 



Let us verify Equation (|39j) for q x . By using (i) as well as the equality x o ip = x, ((39 



is manifestly satisfied on [cj 2 /4, 3w 2 /4] x [0,wi/z[= ^([w 2 /4, 3w 2 /4] x [0, Moreover, 



with (iii) (l39j) is satisfied for q x on ]0,w 2 /4] x [0,ui/i[. Since ip 2 = id, (i39j) is also true for 
q x on [3w 2 /4, w 2 [x [0, and thus finally on the whole C/fi \ [0,u>i[. Likewise, we easily 
verify that (j39j) is valid for q y on C/fi \ ([0, uj\ [+3w 2 /8). Equation (|40|) is immediately 
true, by construction of the continuations. 

It remains to prove that the continuations of q x and (y + l)q y are holomorphic on C/O 
cut along [0,^1 [ and [0, oj\ [+3w 2 /8 respectively. 

We first show that q x is meromorphic on its respective cut parallelogram. The following 
cycles are a priori problematic for q x : [0, uji[, [0, wi[+a; 2 /4 and [0, uj\ [+3cj 2 /4. In an open 
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neighborhood of [0, u)\ [+3^2/4, we have q x (to) = — (y(io)+l)q y (uj)+Q(0, 0, z)+x(uj)y(uj)/z, 
so that q x is in fact meromorphic in the neighborhood of the cycle [0, wi[+3u;2/4. Since 
Equation (f3T)j) holds, q x is also meromorphic near [0,cji[+W2/4 = ip([Q,u)i [+3^2/4) . Thus 
[0, wi[ remains the only singular cycle for q x . Furthermore, q x is clearly holomorphic 
on ]cj 2 /4, 3^2/4] x [0, since it is defined there through the power series Q(x, 0, z). 

On ] 5^2/8, W2[x [0, u>i/i[, we have q x (oj) = —(y(uj) + l)q y (uj) + Q(0,Q,z) + x{u)y{uj) / z 
and the first two terms in the right-hand side of this equality are holomorphic on 
this domain. With Lemma [28] below, the product x{uj)y(uj) may have a pole on this 
domain only at 7^2/8: in fact, x has a pole of order one at this point but y has 
there a zero of order two, so that the product x(cu)y(u}) is holomorphic near 7W2/8. On 
]0, 3w2/8[x [0, uji/z[, we have q x = q x o if), so that q x is holomorphic on this domain, since 
it is on ip(]0, 3w 2 /8[x [0, u>i/i[) =]5w2/8, W2[x [0, it)\/%[. 

A similar reasoning yields that (y + 1)^ is holomorphic on C/f2 \ ([0, oj\ [+3W2/8). □ 

Corollary 27. T/ie function q y is holomorphic on C/O \ ([0, wi[+3w2/8). 

Proof. From Theorem 1251 we know that (y + l)q y is holomorphic on C/O\([0, Wi[+ 3^2/8), 
so that we directly derive that g„ is holomorphic on the same domain, except eventually at 
the points where y + 1 = 0, i.e. at W2/8 and 5^2/8, see Lemma |2"81 However, the generating 
function Q(0, y, z) is bounded at y = —1, see Section^ so that q y (uj) = Q(0, z), being 
meromorphic and bounded near W2/8 and 5W2/8, is actually holomorphic at both these 
points. □ 

The following lemma, which has been used in the proof of Theorem [25| easily follows 
from Lemma[T]and from the fact that on the parallelogram [0, [0, cui/i[, the Weierstrass 
elliptic function p takes each value of C U {00} twice. 

Lemma 28. The only poles of x (of order one) are at CJ2/8, 7^2/8 and its only zeros (of 
order one) are at 3^2/8, 5^2/8. The only pole of y (of order two) is at 3^2/8 and its only 
zero (of order two) is at 7^2/8. The only zeros of y + 1 are at W2/8, 5^2/8. 
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